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1. Introduction 

The renormalised coupling in the Schrodinger functional (SF) scheme has been defined in 
^ and its scale evolution has been studied in QCD with zero and two quark flavours ||3|, We 
here discuss the set-up for studying the running coupling with four quark flavours, in the lattice 
regularisation with staggered quarks [||, |6|]. This writeup is organized as follows. We start with 
a reminder of the basic definition of the renormalised coupling in a formal continuum notation. 
Its lattice formulation with staggered quarks will require a modification at 0(a) of the standard 
framework, even in the pure gauge theory. The corresponding tree-level and one-loop computations 
are then described and we end with an outlook to future work. 

2. The renormalised Schrodinger functional coupling 

The Schrodinger functional is a useful tool to study the scaling properties of QCD. It is defined 
as the Euclidean time evolution kernel for a state at time xq = to another state at Euclidean time 
xq = T. Using the transfer matrix formalism, it can be written as a path integral over fields which 
satisfy Dirichlet boundary conditions at Euclidean times xq = and T and L-periodic boundary 
conditions in space. More precisely, one imposes homogeneous boundary conditions on the quark 
fields, 

(l + 7o)Hvo=o=0=(l-ro)V^Ur 

V^(l-7o)U=o=0='?^(l + ro)Ur' (2-1) 
while the spatial components of the gluon fields satisfy 

^^WL-o=o = Ak{x)\^^^^ = C[. {1.1) 

The SF is considered a functional of these boundary gauge fields, 

^[C',C] = J &[A,\if,v]e'^^^'^'^^- (2.3) 

To define the SF coupling we follow [|l|] and choose Abelian and spatially constant gauge boundary 
fields C and C, 

Ck = Uiag{(^uh,h), C', = Umg{(^[,(^i,(^;). (2.4) 

Under mild conditions on these angular variables, the absolute minimum of the action is uniquely 
determined up to gauge equivalence. Denoting this minimal action configuration by B^{x), one 
may thus unambiguously define the effective action 

r[S] = - In ^[C,C']. (2.5) 

In the weak coupling domain, the SF can be computed by performing a saddle point expansion of 
the functional integral about the induced background field B, leading to an asymptotic series of the 
form, 

T[B] lro[B] +ri [B] +glT2[B] + .... (2.6) 
So 
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The leading term of the series is given by the classical action, while the higher order contributions 
are sums of vacuum bubble Feynman diagrams with an increasing number of loops. 

The finite size scaling technique is based on the idea of a renormalized coupling that runs with 
the box size L. In order for L to be the only external scale in the system all dimensionful parameters 
are taken proportional to L. In particular, one sets T = L and the quark masses are set to zero. Then 
one lets the background field B depend on a dimensionless parameter rj and defines a renormalized 
coupling, 

1 r' 

- ° (2.7) 



g2(L) r" 

where the prime indicates differentiation with respect to t] at Tj = 0, 

r' = ^r[B(r7)]|^^„. (2.8) 

Note that the numerator Fg merely serves as a normalization constant to ensure g = go at tree level. 
When initially defined with the lattice regularisation in place the coupling is well-defined beyond 
perturbation theory. In particular, its non-perturbative running has been previously computed for a 
pure Yang Mills theory in and for QCD with two flavours in [Q]. Here, we prepare the set-up 
for similar studies in QCD with 4 flavours. Staggered fermions seem to be a natural choice, as they 
come in multiples of 4 species, due to the incomplete elimination of the doubler modes. 



3. Subtleties with staggered fermions 

As observed previously in [|5|, ^], the SF with staggered fermions requires lattices where the 
time extent T/a is odd while the spatial extent L/a has to be even. This is illustrated in Fig. [l|. As 
the Dirac spinors are reconstructed from the one-component fields on the corners of a hypercube 
(indicated by different colours) the constraint arises from having the degrees of freedom for a 
multiple of 4 Dirac spinors fit on the lattice, taking into account that half the Dirac spinors satisfy 
Dirichlet conditions at the time boundaries. 

The continuum limit for the SF coupling is usually taken setting T = L already at finite values 
of the lattice spacing. Obviously, with staggered fermions this can only be done up to terms of 
0(fl;). In order to define the continuum limit, it is convenient to imagine a dual lattice, as indicated 
in Fig.[T[ For the spatial directions, the dual lattice has the same number of points as the original 
lattice. However, in the time direction, the number of points is reduced by one. Denoting the 
temporal length of the dual lattice by T' /a one then has T' = L. A further justification may be 
obtained by looking at the fermionic degrees of freedom: the four Dirac spinors are reconstructed 
from the 2^ one-component spinors of the corners of a hypercube and may be assumed to live on 
a lattice with nodes in the center of each hypercube. This is nothing but the dual lattice where 
every second node is eliminated in each direction. In order to get an integer multiple of 4 Dirac 
spinor fields, it is thus necessary that both L/a and T' /a are even. Finally, the case T = L — a can 
be treated the same way, if one imagines a modified dual lattice overlapping by a/2 in both time 
directions. In conclusion, lattices with T = L±a are interpreted as having physical time extent 
T' with T' = T + sa with 5 = ±1, and it is then possible to set T' = L. This modification affects 
even the pure gauge theory, so that we have to revisit the O(fl') effects there before turning to the 
fermionic contributions. 
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Figure 1: 2 fermionic degrees of freedom on a li = 2 dimensional lattice. 



4. Pure gauge theory 

On the lattice we choose the usual Wilson plaquette gauge action, 

S[U] = ^'£w{p)tr{l-U{p)}, (4.1) 
8op 

where the sum runs over all oriented plaquettes p, and U {p) denotes the parallel transporter around 
p. The coefficients w{p) are weight factors to be specified shortly. Due to the abelian nature of the 
boundary fields C and C' the lattice version of Eqs. (2.2) for the link variables U{x,n.) reads 

[/(0,x,A:) =exp(aQ), [/(r,x,A:) = exp(aC[). (4.2) 

With these boundary conditions it is known that the 0(a) effects generated by the presence of the 
boundary are encoded in a single operator, 

d\ii{F^uFok]. (4.3) 

xo=QJ 



An O(a') improved lattice action can thus be obtained from Eq. (11), by setting w{p) equal to 1 
except for the time-like plaquettes attached to one of the time boundaries, where one sets w{P) = 
Ct{go)- The coefficient c, then multiplies a lattice version of (4.3), and, if chosen appropriately, the 
0(a) effects in observables are cancelled. 



4.1 Tree-level considerations 

In perturbation theory, c, is expanded 

cM=cf^+c\'\l + 0{gl). (4.4) 

In the standard framework with T = L one has cf'^ = 1 and the next two coefficients are known, 
too [^. However, when taking the continuum limit at fixed T' /L= 1, a modification of the tree- 
level coefficient cf''^ is required. To calculate it, we first need to determine the minimal action 
configuration as a function of cf^ . 
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4.2 Equations of motion 

In order to be able to write the equations of motion concisely, we follow ^ and introduce the 
covariant divergence of the plaquette, sUghtly modified by the inclusion of the weight factors, 

3 

v=0 

- w[Pf,y{x-v)]uHx-v,v)P^y{x-v)U{x-v,v)}. (4.5) 

The lattice action will be stationary if and only if the traceless antihermitian part of d^P{x,iJ.) 
vanishes, 

d:P{x,^)-d:P^x,^)-^tr {dlP{x,p.)-dlP\x,^)]=0. (4.6) 
We make an ansatz of the form. 



V{x,ix) = exp {aB^{xQ)) 



(4.7) 



with a spatially constant and Abelian field B^{xq). Up to gauge equivalence, the equations of 
motion are then solved by. 



Bo(-^o)=0, Bk{xQ) = < 



Ck xq = 

C( xq = T 



(4.8) 



where / can be computed either numerically or as a power series in a. To check whether the 
above ansatz really leads to the absolute minimum gauge configuration, a cooling procedure has 
been applied to random gauge configurations. While this does not constitute a proof, the apparent 
absence of configurations with lower action for various lattice sizes lends further support to our 
assumptions. 



4.3 Choice of c,^ 



(0) 



The tree-level coefficient c^^^ is to be chosen such that 0(a) effects in observables are can- 
celled. To this order we may just require the lattice action to coincide with its continuum counter- 
part up to 0{a^) terms. Expanding the lattice action to 0(a) 



Sla 



-2 + S + 



Ac) 



(0) 



JO) 



+ 0{a 




(4.9) 



we thus conclude 



JO) 



2 + 5 



(4.10) 



A closer look then reveals that this choice even removes the lattice artefacts in the action up to 
order a^. 
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4.4 One-loop calculation 



Working in a renormalisable gauge, the first two terms in tiie effective action, Eq. (2.5), are 
given by 



ro[B]=gls[V], ri[B] = ^IndetAi -IndetAo 



(4.11) 



where Aq and Ai are the fluctuation operators for the ghost fields and gauge fields respectively. The 
SF coupling to this order then becomes 



f{L)=gl + mi{L/a)gl 



To compute the quantities. 



5^(lndetA,)/r;„ 



mi{L/a) = -r'l/rQ. 



0,1, 



(4.12) 



(4.13) 



we have followed the strategy used in One expects that mi (L/a) has an asymptotic expansion 
of the form, 

oo 

mi (L/a) ''Z-^'' £ (a/L)"(B„ + A„ln(L/a)). 



(4.14) 



n=0 



The results obtained for nii (L) have been confirmed by an independent calculation performed by 
S. Takeda and U. Wolff The (preliminary) results obtained for Bq and Bi are shown in Table |l|, 
where we have set Aq and Ai to their expected values (after having confirmed them numerically). 



s 


Ao 


Bo 


Ai 


Bi 


-1 


22/(471)2 


0.368283(1) 





-0.2318(3) 





22/(471)2 


0.3682817(7) 





-0.1779(3) 


1 


22/(471)2 


0.3682818(7) 





0.1232(4) 



Table 1: Preliminary results for the coefficients in the asymptotic expansion (in collaboration with S. Takeda 
and U. Wolff). 



4.5 Determination of c, 

To determine the one-loop coefficient, we expand the lattice action as a Taylor series about 

(0) 

c c I I 2 d^latt I (1) . . . -V 

Siatt = Siatt]^ JO) +gQ^ — I . '. (4.15) 

Inserting this expansion in the definition of the coupling we arrive at, 

f{L) =gl+ (mi (L/a) - c\'^ d,,T',\^^^^, + ^(^o) (^.16) 

The factor multiplying behaves like 0(a/L), and can thus remove the contribution of Bi, if we 
adjust cP'' properly. We obtain 
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5. Staggered fermion action 



The feimionic part of the action takes the form, 




where the last two terms encode fermionic boundary terms [^]. The coefficient Cf also receives 
a contribution from the fermionic part of the action. We have obtained preliminary results for the 
contribution to mi{L/a); in particular, we find that Aq and Bq coincide with the results obtained by 
Heller [^]. However, before we can quote a value for the staggered one-loop contribution to c|^', a 
more detailed analysis of the fermionic 0(a) boundary counterterm needs to be carried out. 

6. Conclusions 

First steps have been taken in a study of the SF running coupling in QCD with four quark 
flavours. Staggered fermions are a natural choice but require some revision of the standard frame- 
work, due to the fact that L = T can only be imposed up to 0{a) terms. Our proposal to take the 
continuum limit at fixed T' =L with T' = T ±a modifies the O(fl') improvement counterterm pro- 
portional to Cf, which we then determined perturbatively at the tree and the one-loop level. As a 
byproduct this yields an alternative definition of the pure gauge SF which is currently being tested 
in simulations (in collaboration with S. Takeda and U. Wolff). In the near future, this will hopefully 
be followed by numerical simulations including staggered fermions. 
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